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CLASSICAL EIGHTH- AND LOWER-ORDER RUNGE-KUTTA-NY STROM 
FORMULAS WITH A NEW STEPS IZE CONTROL PROCEDURE 
FOR SPECIAL SECOND-ORDER DIFFERENTIAL 
EQUATIONS 


INTRODUCTION 

1. In an earlier report [ 1] this author derived Runge-Kutta-Nystrom (RKN) 
formulas for a special class of second-order (vector) differential equa- 
tions, 

x = f(t,x) , (l) 

which do not contain the first derivative x on the right-hand side. In 
report [ l] , the stepsize control for the Runge-Kutta-Nystrom formulas 
was based on the leading term of the local truncation error in x. 

2. In this report we will derive Runge-Kutta-Nystrom formulas for differen- 
tial equations (l) with stepsize control based on the leading term of the 
local truncation error in x. The derivation of such formulas is similar 
to the derivation in report [ 1] . However, because of the higher accu- 
racy in x required for our new formulas, their derivation will be more 
involved. The new formulas will, in general, require a slightly higher 
number of evaluations of the differential equations per integration step 
compared with the formulas of report [ 1] : for the eighth-order formula, 
13 evaluations instead of 11; for the seventh-order formula, 10 evalua- 
tions instead of 9; and for the sixth-order formula, 8 evaluations 
instead of 7, whereas the fifth-order formula requires 6 evaluations only 
as the corresponding formula of report [ 1]. 

The number of evaluations per step still compares favorably with the 
number of evaluations of the Runge-Kutta formulas for the first-order 
differential equations of report [2]. 

3. The formulas of this report are well suited for problems in which the 
derivative x might assume large values. If x becomes large, the 
stepsize control procedure of report [ l] might lead to somewhat less 
accurate values for x. A stepsize that is small enough for x might, 
for large x, no longer be sufficiently small for x. In Section V an 



example (Orbit 2) is presented in which formulas of this report show a 
higher accuracy than corresponding formulas of report [ 1] . 

Similar to the formulas of report [ 1] , the formulas of this report 
represent pairs of integration formulas for x which differ from one 
another by one additional evaluation of the differential equations. The 
orders of these formulas differ by one. Therefore, the difference of 
the two formulas represents an approximation for the leading term of 
the truncation error in x for the lower-order formula. 

For x, in this report, a Runge-Kutta-Nystrom formula is used, the 
order of which is equal to the order of the higher-order formula for x. 


SECTION I. EIGHTH-ORDER FORMULA RKN 8(9) 


Since the equations of condition for our Runge-Kutta-Nystrom formulas 
are derived in [ l] , the reader is referred to their listing in Table 1 of 
[ 1 ]. 


For an eighth-order Runge-Kutta-Nystrom formula, we allow for 13 
evaluations per step: 


f 0 = f(t 0 ,x 0 ) 


f 

K 



+ a h, x 0 + x 0 a h + 

K K 



x=o 



(k = 1, 2, , 13) 


( 2 ) 


There are actually 14 evaluations. However, as in our earlier report 
[ 1] , we require that the last evaluation f 13 of our differential equations 
can be taken over as first evaluation f 0 for the next step. Therefore, 
only the very first integration step will require 14 evaluations. Our new 
Runge-Kutta-Nystrom formulas are of the form: 


Z c K f K + °( h10 ) 1 

:=0 


x = x 0 + x 0 h + h 2 • 


(3) 



12 


x = x 0 + h • Yj ° K f K + °( h9 ) 
K = 0 


X = x 0 + h * Yl *K f K + °( hl °) 
K = 0 


(3) 

(con. ) 


The quantities t 0 , x 0 , x 0 in (2) and (3) are the initial values for the 
integration step under consideration, while h stands for the integration 

. ❖ 

stepsize. The coefficients a , y ., c , c , and c must now be deter- 
^ k k\ k k k 

mined in such a way that the first equation (3) is a ninth-order approxi- 
mation of x and the last two equations (3) are eighth- and ninth-order 
approximations of x. 


Similar to report [ l], we assume that 

(k = 0, 1, 2, , ll) 


* 

c = c 

K K 


c 12 = 0 


c 13 “ c 12 


(4) 


Because of our assumption for the 14th evaluation (2), the following 
relations must hold: 


y i3A C X 


(A = o, 1, 2, ... , 12) 


a 12 = = 1 


13 


(5) 


We now have to determine the coefficients a , y c , c of formulas 

K K A K K 

(2) and (3). We start with the computation of the weight coefficients 
c and c . From Table 1 of [ 1] we obtain, by assuming that 

K K 


Cl = 

= C 2 = C 3 

= c 4 = 

C 5 

= c 6 = 

0 

Cl : 

= C 2 = c 3 

= c 4 = 

C 5 

— — 

0 


( 6 ) 


3 



the following equations of condition for the coefficients c 
respectively: 

1 

c 7 a 7 + c 8 ff 8 + c 9 £*9 + c 10 af 10 + c u a u + c 12 = - 

c 7 a 7 + C 8 «8 + c 9 a 9 + c 10 a 10 + Cu + c i 2 = ~ 

c 7 a 3 7 + c 8 a 3 g + c 9 a| + c 10 a ? 0 + Cp a?i + c i 2 = ^ 

c 7 a \ + c 8 ffg + c 9 ag + c 10 cvjo + c n Qfu + c l2 = ^ 

c 7 a 7 + Cg a\ + Cg a\ + C 10 a 5 10 + c n a\ t + c 12 = ~ 

c 7 af + c 8 «| + c 9 <a§ + c 10 afi 0 + c^ + c 12 = -r- 

ob 

c 7 ® 7 "t Cg Q?g + Cg 0?g + Cjq Ojq + Cjj Q^u + Cj 2 — 


c 7 a 7 + c 8 ff 8 + eg a 9 + c 10 a 10 + c^ + c 12 


1 

2 


c 7 a 7 + c 8 oi\ + Cg Qf| + c 10 a\ 0 + c n a\ { + c 12 
c 7 ce 3 + c 8 a$ + c 9 a | + c 10 a 3 i0 + c n a? t + c 12 


1. 

3 

j. 

4 


c 7 a\ + c 8 «g + c 9 a\ + c 10 aj 0 + c u c>4i + c 12 - — 


c 7 ® 7 + Cg C?g + Cg Og + Cjq 0!jq + Cj^ ®u + C^ 2 


1 

6 


c 7 af + Cg 0 f| + c 9 a§ + c 10 + c n + c 


12 


1 

7 


c 7 a\ + c 8 ff^+ eg cv^ + c 10 ai 0 + c n a\ x + c 12 = | 
c 7 a? + c 8 a;jj + eg ag + c 10 a \ 0 + c n afj. + c 12 = i 


and c , 

K K 



4 



Equations (7) and (8) must hold if the right-hand sides of (3) are to be 

A 

ninth-order approximations of x or x . 

For given coefficients g , system (7) represents a system of seven 

K 

equations for the six coefficients c 2 , c 8 , c 9 , c 19 , c^, and c 12 . Since there 
is one equation more than there are coefficients c^, there must be a 

restricting condition between the coefficients g^. 


By eliminating from (7) the coefficients c^. Cu, Ciq, • • • one after the 
other, we obtain the following condition: 


Ni 


Dl 


= I . Ni 
3 D t 

= 84 g 8 a 9 a 10 a n - 42 (o; 8 g 9 g 10 + g 8 g 9 a n + g 8 g 10 a n + g 9 

+ 24 (a 8 g 9 + g 8 g 10 + « 8 a ll + a d “10 + “11 + «lo «ii) 

- 15 (g 8 + Gg + CtjQ + Gjj) + 10 

= 70 a 8 a 9 a 10 a n - 28 (g 8 g 9 g 10 + a 8 a 9 g u + g 8 g 10 + g 9 

+ 14 (g 8 a 9 + a 8 g 10 + a 8 G n + a 9 a 10 + a 9 a lt + G 10 G tl ) 

_ 8 (g 8 + a 9 + a 10 + %) + 5 


“to 


“to 



( 9 ) 


Similarly, system (8) yields the following two restricting conditions: 


n 2 


+ 


1 N 2 

2 ' D 2 

70 g 8 a 9 a 10 g u - 42 (a 8 g 9 a 10 + a 8 a 9 a n + g 8 g u + g 9 a 10 a n ) 


28 (a 8 a 9 + a 8 a 10 + a 8 a n + g 9 a 10 + g 9 a n + g 10 g u ) 
20 (g 8 + g 9 + g 10 + g u ) + 15 


(10) 


5 



D 2 = 70 <* 8 a 9 a 10 a n - 35 (a 8 a 9 a 10 + a 8 a 9 a n + a 9 a 10 a n + a 9 a 10 Qf u 
+ 21 (« 8 oi 9 + a 8 a 10 + a 8 a n + a 9 a 10 + a 9 a n + a u a lt ) 

- 14 (Q! 8 + a 9 + a 10 + a lt ) + 10 

and 



N 3 = 126 a 8 « 9 a 10 «u - 84 (a 8 a 9 a i0 + a 8 a 9 a n + a 8 a 10 a n + a 9 a 10 a 

+ 60 (or 8 a 9 + a 8 a 10 + o? 8 a n + a 9 a 10 + a 9 a n + a u o^) 

_ 45 (a 8 + o; 9 + a 10 + Qf u ) + 35 
D 3 = N 2 

Assuming that 

+ « 8 = 1 , a io + a n = 1 » 

the conditions (9), (10), and (ll) reduce to 

i 84 o? 8 a 9 of 10 a n - 18 (a 8 a 9 + a 10 a tl ) + 4 

7 ~ 3 70 a 8 a 9 a 10 a u - 14 (a 8 a 9 + a 10 ar u ) + 3 

a* = — 

1 2 

x 126 a 8 a 9 a 10 a n _ 24 (ar 8 a 9 + a 10 a lt ) + 5* 

7 3 70 a 8 a 9 ar 10 a lt - 14 (a 8 a 9 + ar 10 a u ) + 3 

6 


I 


) 

( 10 ) 

(con.) 


) 

( 11 ) 


( 12 ) 

(13) 



Inserting ce 7 = — into the first and third equation ( 13) yields both times 
2 

the following relation between o; 10 Q^i &od ^8 a 9 : 


«io a n 


l . 6 o; a Q! 9 - 1 
6 7 a 8 o!9 - 1 


(14) 


1 

We select o' 8 = — , 

O 


9 

— as two values compatible with the first equa- 
3 


tion (12) and leading to relatively reasonable values for the remaining 
coefficients of our formula. 


From (14) we then obtain a i0 a n = and because of the second equa- 
tion (12), = (5 - nTi 5 ) , Q; ii = Yo ( 5 + ,N ^ - ^)* 


Listing all the coefficients a that enter equations (7) and (8), we have 

K 


a l = ^ , a 8~3> “9 - f » “iO - Yo ( 5 " ) * 

( 5 + n/~15 ) . 


(15) 


For these coefficients a , equations (7) and (8) can be solved with 

K 

respect to c and c • We find that 
r K K 


16 _ 243 _ 243 

C? " 315 ’ ° 8 1540 ’ C& 3080 ’ 10 

Cli = 1386 » C 12 = 0 


25 

1386 


(5 + n/ 15 ) , 


( 16 ) 


and 


7 



by 


32 

315 ’ CfS 


729 . = 729 . = 125 

3080 ’ ° 9 3080 ’ Cl ° 693 


125 . _ _9_ 

c ii 693 » c 12 280 


( 17 ) 


7. Next we have to find the coefficients y .. Similarly as in [1], we 

K A 

assume for the y . the following relations: 

kX 


K-l 


Z r - a * = “ a 


1 3 

= — rv a 


A=1 


r «A A 6 K 


(k = 2, 3, 4, , 12) 


K-l 

Z 

A=1 


2 -L 4 

v cr = — or 
icA A 12 K 


(k = 2, 3, 4, ... , 12) 


K-l 

2 

A=1 


3 1 5 

y cr = — cr 
Va A 20 k 


(k = 4, 5, 6, ... , 12) 


K-l 

Z 

A^l 


y oA. = — oA 
r K X X 30 K 


(k = 5, 6, 7, ... , 12) 


We further assume that 

751 = Tei = Yn = Ysi = Y91 = Yioi = Tin = Y121 = 0 

^72 “ ^82 = ^92 = Yl02 “ Yil2 = Yl22 = 0 

773 = 783 = ^93 = Yl03 = Yll3 = Yl23 = ® 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


8 


l 



and 


11 

Tj c 7=0 
u „ K k4 
k=7 


11 


E ®.. t., + c -,o \ 


K=1 


K k4 12 


124 


> = 0 


(23) 


The upper line in (23) holds for the eighth-order formula and the lower 
line for the ninth-order formula. Because of c 4 = 0 , we obtain from 
(23) 


Tm = 0 


(24) 


As pointed out in [ l] , such assumptions as (18) to (2l) and (23) convert 
the necessary and sufficient equations of condition (Table 1 of [ l] ) for 
the coefficients y into a system of sufficient equations of condition that 

K A 

can be solved in a relatively easy manner. 


If we proceed in a way quite similar to that of [ 1] , the application of 
(6) and (18) to (23) reduces Table 1 of [1] to 


1 

c 0 + c 7 + c 8 + c 9 + c 10 + c 14 - — 
c 0 + c 7 + Cg + c 9 + c 10 + c lt + C 12 = 1 


(25) 


to (7), (8), and to the following four equations: 


c 7 Pft + c 8 Pg 5 + c 9 P 95 + c 10 P 106 + c n P 115 - 3024 


^7 -^75 + ^8 ^85 + ^9 ^*95 + ^10 ^*105 + ^11 ^*115 + ^12 



(26) 


1 

336 


9 



c 7 c * 7 P 75 + c 8 a 8 P 85 + c 9 a 9 P 95 + c 10 a 10 P 105 


+ cn «il Pu 5 + c 12 


_1_ 

42 


1 

378 


°7 P 76 + ^8 P 86 + P 96 + ^10 P 106 + ^11 P 116 + ^12 ^ ~ 


(26) 
(con. ) 


In (26) we have used the abbreviations 


v . v , 

y a + y a + 

>1 1 >2 2 


+ y . oi = p 

nn-l n-l nv 


(27) 


Equations (26) correspond to equations (IX, 27), (VIII,27) , (IX, 43)*, 
and (IX, 53)’ of Table 1 in [ 1] . 

From the second equation (26), we obtain 


p i25 " 42 


(28) 


The coefficients y can now be determined from equations (18), (19), 

K A. 

(20), (21), (23), and (26). 

8. There are restrictive conditions for the a f s that we derive in the 

K 

following. From equations (18) and (19) _ , we obtain 

K 2 K — 2 


a i = l a 2 • (29) 

From equations (18) , (19) , (20) , and (21) , the following 

K O K 0 K 0 K — 5 

condition results: 


10 


i 



( 30 ) 


5 a 3 a 4 ~ 3 (a s + a 4 ) a 5 + 2 a\ 

a 2 = a 5 

10 a 3 a 4 - 5 (a 3 + a 4 ) a 5 + 3 a| 

and from equations (18) k=? , (19) k=? , (20)^, and (21)^, 

5 a 5 a g — 3(a 5 + a g ) a 7 + 2 oi\ 

a 4 = a 7 • . ( 31 ) 

10 a s a 6 - 5(a 5 + ag) a 7 + 3 a| 

Having selected a 7 by equation ( 15 ) , the coefficients a 3 , a 5 , and ag 
can be chosen arbitrarily, while the coefficients a 4 , a 2 , and a t are 
determined by equations (3l), (30), and (29). 

9. Since the expressions for the coefficients y . can be derived in a 

kX 

straightforward manner, we only sketch their derivation: 
from (18) , we obtain y 21 ; 

K Z 

from (18) k=3 , (19) k=3 , we obtain y 31 , y 32 ; 

from (18) k=4 , (19)^, (20) k=4 , we obtain y 41 , y 42 , y 43 ; 

from (l8) K=g , (19) k=5 , ( 2 °) k=5 » we obtain y 52 , y 53 , y 54 ; 

from (18) k=6 , (19) k=6 » ( 2 °) k=6 * ( 21 ) k=6 > we obtain y 62 , y 63 , yg 4 , y 65 ; 

from (18) k=7 , (19) k=? , ( 2 °) k=7 » we obtain y 74 , y^, y 76 ; 


from (l8) K=g , ( 19)^ =g , (20) K=g , ( 21 ) (f=8 » we obtain y M , ygg, y M , y 87 . 

With 

T94 = 0 , (32) 


11 



we find from 


( 18 )k= 9 ’ ( 19 ) k=9 ’ ( 2 °) k = 9 ’ ( 21 ) k=9 the coefficients y 95 , y 96 , y 97 , y 98 . 

Now we can consider (23) as a system of two equations for the two 
coefficients y 104 and y 114 . 

The first three equations (26) can serve to express P 105 and P 115 by P^, 
p 85» p 95 and the weight factors c and c . Using the values (16) and (17) 

for the weight factors, we find that these first three equations (26) are 
not independent: The second equation is equal to the sum of the first and 
the third equation. Therefore, we may omit one of these three equations 
and obtain P 10B and P 115 from the remaining two equations. 

Since Y m and P 105 are now known, equations (18) , ( 19) _ , 

( 2 °) k== 10» ( 21 ) k =io’ and the e 9 uation ( 2? ) defining P 105 yield y 105 , y 106 , 

Tl07» 7l08> an d 7109* 

From the fourth equation (26), we obtain P U6 as a function of P 76 , P 86 , 
p 96> p i06 and the weight factors c . The values for y 114 , P 115 , and P 116 

being known, we find y 115 , y 116 , y 117 , y 118 , y 119 , and y 1110 from equations 
( 18 ) k=11 » ( 19 ) k=11 > ( 2 °) k=11 ’ ( 2 1) k=11 and equation (27) for P 115 

and P 116 . 

Putting 


7125 = 0 ( 33 ) 

and using (28), we have five equations: (18) _ , (19) , (20) 

K- 12 k=12 ' k =12 

( 21 ^ = 12 ’ and the definition (27) for P 125 to express y 127 , y 128 , y 129 , y 1210 , 
and y 1211 by y 126 . 

10. The coefficient y l26 remains undetermined. Collecting from Table 1 of 
[ lj all terms that contribute to the leading term of the truncation error 
for the second equation (3), we find that for given coefficients a this 

K 

leading term is proportional to y 12e . Therefore, by varying y 12e we can 


12 


1 



influence the local truncation error in x. However, we must not choose 
Ym = 0, because our second equation (3) would then become a ninth- 
order formula and our stepsize control procedure would break down. 

11. The coefficients y are determined from the well-known relation, 

K 0 


K-l 


E 

AfO 


y 


K A 


1 

2 


ct 


J 


(34) 


and the weight factors c 0 and c 0 from (25). 


12. This terminates the computation of the coefficients of our eighth-order 
Runge-Kutta-Ny strom formula. We shall label our formula RKN 8(9)- 
13-x to indicate that we have a pair of formulas of eighth- and ninth-order 
(in x) with 13 evaluations per step and with stepsize control based on the 
leading term of the local truncation error in x. 


In Table 1 the pattern of such a formula RKN 8(9)-13-x is presented. 
All coefficients different from 0 and from 1 are indicated by asterisks 
(*). Comparing Table 1 with Tables 5, 8, and 11 of the lower-order 
formulas makes it easy to recognize the structure of our formulas. 


In Table 2 are listed the coefficients of a formula RKN 8(9)-13-x in 
fraction form. The parameters a 3 , a 5 , and Q! 6 were chosen as 1, and 
— , respectively, since this choice leads to coefficients that are still 
relatively reasonable. 


In Table 3 are listed the coefficients of Table 2 in decimal form (32 

5 

decimal digits). In Table 3 we set = since this value yielded 

the most accurate results for the examples that we ran on the computer. 
Table 4 represents a flow chart for programming our formula 
RKN 8(9) - 13-x. 


13 



SECTION II. SEVENTH-ORDER FORMULA RKN 7(8) 


13. Since seventh-order Runge-Kutta-Ny strom formulas are derived in a 
very similar way as the eighth-order formula of Section I, we restrict 
ourselves to briefly outlining their derivation without presenting all 
details. 


We allow 10 evaluations per step for the seventh-order formula. 
Assuming 


«9 = 

«10 = 1 




Cl = 

c 2 = c 3 = 

: c 4 

= c 5 

= 0 

Cl = 

C2 = C3 = 

: c 4 

= C S 

= 0 


(35) 


we obtain for the weight factors c and c the following equations of 
condition: 


c 6 a 6 + c 7 a 7 + c 8 + c 9 = - 


c 6 a\ + c 7 a* + c 8 a\ + c 9 = — 
c 6 a 3 6 + c 7 a\ + c 8 + c 9 = ^ 


c 6 a\ + c 7 aj + c 8 a\ + c 9 = ^ 


(36) 


c 6 a\ + c 7 a] + c 8 a\ + c 9 = -^ 

Cg a\ + C 7 CXy + Cg Q!g + Cg = ^ 


14 


I 



and 


c 6 a e + c 7 a 7 + c 8 a 8 + c 9 = - 


c 6 a\ + c 7 <4 + c 8 a\ + c 9 


c 6 a z G + c 7 a] + c 8 al + c 9 


c 6 a\ + c 7 a* + c 8 a\ + c 9 


c 6 a\ + c 7 ot\ + c 8 oi\ + c 9 


c 6 aj + c 7 a? + c 8 al + c 9 = - 


c 6 c*l + c 7 a] + c 8 a\ + c 9 


Equations (36) lead to the following restrictive conditions for the cn ’s: 

K 


a e = 


- 14 a 7 a; 8 - 7 (a 1 + o; 8 ) + 4 

7 

5 a 7 o; 8 _ 2 (a 7 + a; 8 ) + 1 


a R ~ — 

6 2 


. 14 a 7 a 8 - 8 (a 7 + a 8 ) + 5 


14 a 7 a 8 - 7 (a 7 + a 8 ) + 4 


whereas from (37) the conditions 


<*6 = 


5 a 7 a 8 - 3 (a 7 + a 8 ) + 2 
10 ct 7 a 8 _ 5 (a 7 + Qf 8 ) + 3 


a a - — 
6 7 


21 a 7 a 8 - 14 (a 7 + Q! 8 ) + 10 


5 a 7 a 8 - 3 (a 7 + o; 8 ) + 2 


15 



28 a 7 o? 8 - 20 (a 7 + a 8 ) + 15 


(39) 
(con. ) 


a 6 ~ 


1 

2 


21 a 7 o? 8 - 14 (a 7 + ctg) + 10 


result. 

Assuming that 

<a 7 + a 8 = 1 , 

equations (38) and (39) reduce to 

a l «8 = | 

and 

1 

a~ = — 

6 2 

From (40) and (42) we find that 

a 7 = (7 -\T21 ) , q-8 = ^ (7 +'/YT) . 


(40) 


(41) 


(42) 


(43) 


The a -values (42) and (43) lead to the following values for the weight 

K 

factors 


c, = — , c 7 = — (7 + \Z~2l) , Co = — (7 - V~2l) , 

6 45 ’ 7 360 V 7,8 360 V 7 

c 9 = 0 


(44) 




16 

45 


c 7 


49 . _ _49_ . _ J_ 

180 ’ Cs ” 180 ’ C9 _ 20 ' 


(45) 


16 


I 



14. Similar to our assumptions in Section I, we now assume that the following 

relations for the coefficients y . hold: 

K X 


K - 1 

Z 

\=1 


-v q; = _ Qf ' 5 

KX X 6 K 


(k = 2, 3, 4 9) , (46) 


K-l 

z 

K = 1 


y Of 
/c\ A 


2 _ 


1 4 

CT 

12 K 


(k = 2, 3, 4, . . . , 9) , (47) 


K - 1 

z 

K = 1 


(«=4. 5, 6 9) 


(48) 


We further assume that 


741 75 1 Tsi _ 7?i ~ 781 “ 791 ~ 0 


702 _ 772 “ 782 ~ 792 ~ 0 


and 


(49) 


c 6 763 + ^7 y 73 + C 8 783 “ 0 


(50) 


These assumptions reduce the equations of condition for 7 listed in 
Table 1 of [ 1] to the following four equations: K 


c 6 P 64 + c 7 P 74 + c 8 P 84 = 


1680 


94 


P 64 + C 7 P74 + C 8 Pg 4 + C 9 J l 

90 


1 

210 


(51) 
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C 6 a s P 64 + C 7 « 7 P 74 + C 8 a 8 Pg4 + C 9 * ^ ^ 

P 65 + ^7 P 75 + ^8 P 85 + ^9 * 42 = 33 ^ 


(51) 
(con. ) 


These equations correspond to equations (VIII, 15), (VII, 15)’, (VM, 
21)’, and (VIII, 27)’ of Table 1 in [ 1] . 

From the second equation (51), it follows that 


P 94 


90 


(52) 


15. 


Again, there are restrictive conditions for the a 's. From equations 
( 46 >/c= 2 and ^ 47 ^k= 2’ WG obtain 


a 


l 



9 


and from equations (46) - (47) A , and (48) 

K= 4 k= 4 k= 4 


(53) 


a 


2 


1 5»q - 3 <a 4 

5 2a 3 - a 4 


(54) 


We also assume that 


“5 = “8 


(55) 


The coefficients a 3 and o; d can be chosen arbitrarily while the coefficients 
a 2 , a lt and o; 5 are determined by (54), (53), and (55). 
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16. We now compute the coefficients y .. The coefficients y 21 , 73i> 732> 

K A 

Y42> 743, 752 > 753 » 754 . 7 6 3 > 764» and 7 6 s are computed from equations (46), 
(47), and (48) in the same way as they were computed in Section I from 
(18), (19), (20), and (21). The only difference is that y 41 and y 62 are 
now zero. 

Inserting (44) and (45) into the first three equations (51), we find that 
these three equations are not independent of each other: The second 
equation (51) is equal to the sum of the first and the third equation 
(51). Therefore, we can omit one of the three equations (5l) and use 
the other two equations (51) to express P 74 and P&J by P^. 

From (46) k=? , (47) K=r?r , (48) k=7 , P 74 , we obtain y 73 , y 74 , y^, y 76 . 

The coefficient Y33 is determined by (50). 

From (46) k=8 , (47) K=sg , (48)^, P M , we obtain y M , y B , y 86 , y 87 . 

By the special choice (55) of a 5 , the fourth equation (51) is also 
satisfied. 

Setting 


795 = 0 , (56) 

equations 

( 4 6)^_ 9 , (47) , (48) , P94, yield y M> 7 96 . 7 97 » 7 98 as functions of y 93 . 

The leading term of the local truncation error in our seventh-order 
formula for x is proportional to y 93 . Therefore, by varying y 93 , we can 
again influence the local truncation error in x. 

17. Tables 5, 6, and 7 show the pattern for a formula RKN 7(8)-10-x and 

the coefficients for such a formula in fraction and in decimal form. For 

1 5 

the free parameter o; 3 and a 4 , we chose — and — since we then obtain 
relatively simple values for the remaining coefficients of our formula. 

By setting 7 93 = jr-, we got the best accuracy in the examples that we 
ran on the computer. 
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For the seventh- and lower-order formulas we do not include flow charts 
because they would be very similar to the flow chart for the eighth-order 
formula (Table 4). 


SECTION III. SIXTH-ORDER FORMULA RKN 6(7) 


18. Again, we only briefly outline the derivation of our sixth-order formula, 
which requires eight evaluations per step. Assuming that 


C 1 " C 2 - C 3 - 0 
Cl = c 2 = C 3 = 0 



(57) 


the weight factors c and c are obtained from 

K K 


i 

c 4 a 4 + c 5 “5 + c 6 + c 7 “ ^ 

c 4 a\ + c 5 al + c s a\ + c 7 = ^ 

C 4 a 4 3 + c 5 a 3 + c s a 3 s + C 7 = ~ 

c 4 oif + c R of + c 6 Qfg + c 7 = — 

c 4 « 4 + c 5 ffj+Cj a^ + c, = j- 


(58) 


and 

c 4 a A + c 5 % + c 6 a 6 + e 7 
c 4 a\ + c s + c s a| + c 7 
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(59) 
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c 4 Ol\ + c s 0's + C 6 0 i\ + 6 7 = i 

c 4 + C 5 a| + c 6 aj + c 7 = | 

c 4 a\ + c 5 a| + c 6 a\ + 6 7 = 7 

c 4 af + c 5 +6 6 aj? + c 7 = | 

Equations (58) lead to the restrictive condition, 

14 a; 5 q; 6 _ 7 (a s + ttg) + 4 

a 4 = - , 

5 » 5 ffg - 2 ( cr 5 + Q!g) + 1 

and equations (59) to the two restrictive conditions, 

5 tVg ctg - 3 (o; s + ttg) + 2 
«4 = 

10 0! s Og _ 5 (a 5 + ag) + 3 

1 21 e* 5 a s - 14 (<* 5 + a 6 ) + 10 
“4 = 7 

5 “s a 8 - 3 ("5 + a 6) + 2 

Assuming that 

«S + a 6 = 1 

equations (60) and (6l) simplify to 




(59) 
(con. ) 
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1 _ 14 a 5 o: fi - 3 
7 5 a 5 ctg _ 1 

1 


(63) 


1 ' 21 o? 5 - 4 

7 5 a 5 a 6 _ 1 


From (63), we obtain 


«S 


a 6 


1 

7 


(64) 


Equations (62), (63), and (64) yield 

" 4 = 2 ’ = lj( 7 <* 6 = -j£ (7+^lT) . (65) 


Equations (58) and (59) then lead to 

c ‘ = 45 • ° 5 = 35o ( 7 + n '^ )> c « = 360 (7 -' rn)l c ’ = 0 ■ 

. 16 . 49 .49 . 1 

° 4 45 ’ ° 5 180 ’ ° 6 180 ’ ° 7 20 


( 66 ) 

(67) 


19. 


Again, we assume that certain relations hold for the coefficients 



K-l 

Z 

A=1 


K A A 


= — or 
6 K 


(k = 2, 3, . . . , 7) , 


( 68 ) 
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(69) 


K - 1 

z 

AfI 


\x “1 = 


1 4 

— or 

12 K 


(k 2, 3, . . . , 7) , 


K-l 

z 

AfI 


3 1 5 

T' cr = — or 
r K\ X 20 K 


(k = 4, 5, 6, 7) 


(70) 


We further assume that 


741 = 7 5 i = 7 6 i = 7 7 i=0 


(71) 


If assumptions (68) to (7l) hold. Table 1 of [ 1] reduces to one equation 

for y ^ : 

K A 


C 4 P 44 + C 5 P 54 + C 6 P 64 + C 7 


30 


1 

210 


(72) 


corresponding to equation (VII, 15)* of this table. 

20. From (68) and (69) we obtain the restrictive condition 

K = 2 K= 2 


a 


l 


1 

2 


a 


2 


and from (68) ., (69) ., (70) , 

v / «=4 ' k =4 k = 4 


(73) 


Of 


2 


1 5 Q!q - 3 

— O', “ 

5 4 2 a 3 - 0^4 


(74) 


while the coefficient a 3 can be chosen arbitrarily. 
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21. The computation of the coefficients y follows the same pattern as in 

Sections I and II. The values for y 21 , y 31 , y 32 , y 42 » Y43* Y52» Ys3» and Y54 
are computed from (68), (69), and (70). For the computation of y 62 , 
y 63 , ygj, y 65 we also use Pg 4 , as computed from (72). 

Putting 


y 72 ~ ° 


(75) 


we can express y 74 , y 75 , y 76 by y 73 using equations (68), (69), and (70) 
with k=7. 

The leading term of the local truncation error in our sixth-order formula 
for x is proportional to y 73 . 

22. In Tables 8, 9, and 10, we present the pattern and the numerical values 
for the coefficients of a sixth-order formula RKN 6(7) -8-x. For our 

formula, we chose a 3 = ^ and y 73 = . 

SECTION IV. FIFTH-ORDER FORMULA RKN 5(6) 


23. We allow for six evaluations per step. Assuming that 


QL — Qfg — 1 

Ci = c 2 = 0 
Ci = c 2 = 0 


(76) 


we find the following equations of condition for the weight factors c 
and c : K 

K 
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c 3 + c 4 a 4 + c 5 = | 

c 3 a\ + c 4 aj + c 5 = 

c 3 a l + C 4 «4 + c 5 = ^ 

c 3 a i + c 4 at + c 5 = ^ 


( 77 ) 


and 


c 3 a 3 + C 4 a 4 + C 5 = 

c 3 Q?1 + c 4 aj + c 5 = | 

C 3 “3 + c 4 aj + C 5 = | 

c 3 af + c 4 af + c 5 = -i 

o 

c 3 orf + c 4 orf + c 5 = | 


(78) 


Equations (77) lead to the restrictive condition 


a 


3 


2 o? 4 - 1 
5a ( - 2 


(79) 


and equations (78) to the conditions 


«3 


a 


3 


1 . 5 a 4 - 3 
5 2 a 4 _ 1 

3 a 4 - 2 
5 - 3 


(80) 
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Equations (79) and (80) can be written in the form, 


5 dig oi^ - 2 (o; 3 + <a 4 ) + 1-0 
10 a? 3 a?4 _ 5 (a 3 + a 4 ) + 3 = 0 
5 a 3 a i - 3 (a 3 + a 4 ) + 2 = 0 

Assuming that 

a 3 + a 4 = 1 , 
equations (81) reduce to 
a 3 «4 = | 

Equations (82) and (83) result in 


( 81 ) 


(82) 


(83) 


a 


3 


_ 1 _ 

10 


(5 - 'v/T ) 


Oi 4 


1 

10 


(5 +\T 5") 


(84) 


For the weight factors c and c , we find from (77) and (78) that 

K K 


c 3 - TT ( 5 + ^ )» C 4 “ ^7 ( 5 - ’Js)* c s - 0 


. _ 5 • _ 5 . _ 1 

Cs " 12 ’ C4 ” 12 * ° 5 " 12 * 


(85) 


( 86 ) 


24. The following relations are assumed to hold for the coefficients y . : 


26 


1 



t y KX \ = H <“ = 2 - 3 - 4 - 5 > • 


K — 1 


a 2 
X X 


1 4 

— a 4 
12 K 


(k = 3, 4, 5) 


c 3 731 + c 4 7 4 i = 0 . 


Table 1 of [1] reduces again to one equation for y 


Cj P 3S + Cl p 43 + h ■ ^ = -go • 


corresponding to equation (VI, 8) of this table. 

From (87) _ , we obtain y 2i as a function of 04 and ol 2 . Similarly 

K 2 

(87) _ and (88) yield y 31 and y 32 as functions of 04, a 2 , and Qf 3 . 
The four equations (87) , (88)^, (89), and (90) for the three 

coefficients y 41 , y 42 , y 43 lead to a restrictive condition for the s: 


a 2 = 


«3 ( a 3 + “l) " 


2 ot z a\ + a\ - a\ 


Inserting (84) into (91) leads to 


a 2 = Yo (5 + ) - 04 


Putting 


752 = 0 
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we use equations (87) _ and (88) __ to express y 53 and y 54 by a u a z , 

K— b K— b 

a 4 , and y 51 . 

Since the leading term of the local truncation error in x is proportional 
to y 51 , we can achieve a small truncation error by a proper choice of 

Ysi- 


25. Tables 11, 12, and 13 show the pattern and the numerical values for 
the coefficients of a fifth-order formula RKN 5(6)-6-x. For our 
formula, we chose = 1 and y 51 = 1/300. 

SECTION V. AN APPLICATION: NUMERICAL INTEGRATION 

OF CENTRAL ORB ITS 


26. Let us consider a particle moving in the (x, y) -plane under the influ- 
ence of an attracting mass, this mass being located in the origin of the 
(x, y)-coordinate system. With proper units, the motion of the 
particle in the (x, y) -coordinate system is described by the differential 
equations 


x = 


x 

-7 



(r = \ r x 2 + y 2 ) . 


(94) 


If we consider an elliptical orbit, equations (94) are solved by 


x = e + cos ip 
y = b • sin ip 
t = ip + e • sin ip 


(95) 


the last equation (95) being Kepler’s equation. 
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In (95) b stands for the minor semi-axis of the elliptical orbit, e for its 

eccentricity (e = \T 1 - b 2 , the major semi-axis a being equal to l), 
and ip for the eccentric anomaly. 

From (95) we can obtain the coordinates x, y for the particle and the 
time t as functions of ip. Differentiating (95) with respect to t, we 
can also express the velocities x and y by ip: 


x 


sin ip 

1 , V 

1 + e cos ip 


b cos ip 
1 + e cos ip 


(96) 


For ip=0 and ip=2ir (one complete orbit), we find from (95) and (96): 


t = 0 

> x = 

t = 2 tt 


1 + e 


y = 0 , x 


0 , y 


b 

1 + e 


(97) 


27. Integrating (94) for a complete orbit, we can check the accuracy of the 
numerical integration by comparing the results of the integration with 
(97). We applied our RKN formulas to two elliptical orbits, one of 
them (Orbit No. 2) being highly elliptical. The orbits are pictured in 
Figure 1 (page 50). 

For comparison, the same orbits were also computed using the Runge- 
Kutta-Nystrom formulas of report [l] and the Runge-Kutta formulas of 
report [ 2] . 

24 7 

For Orbit No. 1 we chose e = — , leading to b = — and the initial 

25 25 

values 


x 0 = — = 1. 96 ; yo = 0 , X 0 = 0 , y 0 


l 

7 


(98) 


For Orbit No. 2 we chose 


924 

925 


leading to 



and 
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1849 

925 


_ 43 

1.99 891 . . . , y q = 0 , Xg = 0 , Yq ~ ^g^g 


. ( 99 ) 


Xu — 


In Tables 14 and 15 the results of the various integration formulas are 
listed for Orbits Nos. 1 and 2. The tables show for a complete revolu- 
tion (t=27r) the number of integration steps, the errors, and the execu- 
tion time of the various methods. All computations were performed on 
an IBM-7094 computer in double precision (16 decimal digits); all 
computer programs were written in IBM machine language. 

For all formulas, the stepsize control was applied as described in No. 

26 of report [ 1] , using a tolerance of 0. 1 • 10 -16 . 

28. Tables 14 and 15 show clearly that, for these examples, the Runge- 
Kutta-Nystrom formulas of this report and those of report [1] are 
considerably faster than our Runge-Kutta formulas of report [2]. The 
fifth-order formula of this report is also faster than the fifth-order 
formula of [1]. For higher orders both types of Runge-Kutta- Nystrom 
formulas required, in these examples, approximately the same time on 
the computer. However, for Orbit No. 2, the seventh- and eighth-order 
formulas of this report yield more accurate results than the correspond- 
ing formulas of [ 1] . It is believed that this increase in accuracy is 
gained by the new procedure of stepsize control that we have applied in 
this report. In the vicinity of singularities it seems to be more reason- 
able to base the stepsize control on the x- values than on the x- values. 
Orbit No. 2 comes rather close to the attracting mass (x = y = o), and 
the particle moves fast in the vicinity of the mass. Therefore, in this 
case, a stepsize control based on the x-values seems to be more 
appropriate to the problem. 


Computation Laboratory 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812, February 25, 1973 
014-00-00-00-62 
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TABLE 2. COEFFICIENTS (IN FRACTION FORM) FOR RKN 8(9)-13-x 
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TABLE 3. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 8(9)-13-x 


“l 

= 

0. 3333 

3333 

3333 

3333 

3333 

3333 

3333 

3333 


“2 

= 

0. 6666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 


“3 

= 

0.5 









“4 

= 

0. 3333 

3333 

3333 

3333 

3333 

3333 

3333 

3333 


“s 

= 

1 









“6 

= 

0. 1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 


a i 

= 

0. 5 









“8 

= 

0. 3333 

3333 

3333 

3333 

3333 

3333 

3333 

3333 


“9 

= 

0. 6666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 


“10 

= 

0. 1127 

0166 

5379 

2583 

1148 

2073 

4600 

2176 


“11 

= 

0.8872 

9833 

4620 

7416 

8851 

7926 

5399 

7824 


“l2 

= 

1 









“13 

= 

1 









Tio 

= 

0.5555 

5555 

5555 

5555 

5555 

5555 

5555 

5556 

10- 1 

T 20 

= 

0. 7407 

4074 

0740 

7407 

4074 

0740 

7407 

4074 

10" 1 

7n 

= 

0. 1481 

4814 

8148 

1481 

4814 

8148 

1481 

4815 


T 30 

= 

0. 5468 

75 







10- 1 

Y*i 

= 

0.7812 

5 







10 -1 

Y 32 

=- 

-0. 7812 

5 







10" 2 

740 

= 

0. 2746 

9135 

8024 

6913 

5802 

4691 

3580 

2469 

10" 1 

7 4 l 

= 

0. 5740 

7407 

4074 

0740 

7407 

4074 

0740 

7407 

10 -i 

742 

= 

0. 1018 

5185 

1851 

8518 

5185 

1851 

8518 

5185 

10- 1 

743 

=- 

■ 0. 3950 

6172 

8395 

0617 

2839 

5061 

7283 

9506 

10- 1 

7 5 o 

= 

0.9166 

6666 

6666 

6666 

6666 

6666 

6666 

6667 

10 _1 

7 5 1 

= 

0 









752 

= 

0.225 









753 


•0. 2666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 


7 5 4 

= 

0.45 









7 6 o 

= 

0.4693 

7007 

1183 

8190 

8810 

1971 

8087 

2383 • 

10“ 2 

7 6 i 

= 

0 









7 S2 

= 

0.2178 

2756 

6935 

9345 

6282 

0708 

2253 

7215 * 

10 -2 

763 

==- 

■ 0.5315 

3595 

6013 

9319 

3223 

7068 

6492 

0095 • 

10 -2 

764 

= 

0.4773 

1544 

9880 

6076 

3095 

0566 

4786 

8719 • 

10 -2 

7 6 5 


■0. 1569 

3181 

3691 

4539 

9018 

8939 

1296 

4939 • 

10" 3 

7 ? 0 

= 

0. 1510 

4166 

6666 

6666 

6666 

6666 

6666 

6667 • 

10 _1 

7?l 

= 

0 









772 

= 

0 









773 

= 

0 








10 _1 

774 

= 

0.3867 

1875 






• 

775 

= 

0.3255 

2083 

3333 

3333 

3333 

3333 

3333 

3333 * 

10 -4 


34 


I 




TABLE 3. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 8(9)-13-x 

(Continued) 

y 76 = 0.7119 1406 25 • 10" 1 

y 80 = 0. 1069 9588 4773 6625 5144 0329 2181 0700 • 10" 1 
?8i = 0 

782 = 0 

783 = 0 

784 = 0.6944 4444 4444 4444 4444 4444 4444 4444 • 10~ 2 

y® = 0. 1286 0082 3045 2674 8971 1934 1563 7860 • 10“ 4 

ygg = 0. 3883 9285 7142 8571 4285 7142 8571 4286 * 10" 1 

y 87 =-0.9406 2316 2845 3850 6760 7289 8295 1205 • 10“ 3 

y 90 = 0. 1810 6995 8847 7366 2551 4403 2921 8107 • 10 _1 

7 9 1 = 0 

7 92 = 0 

793 = 0 

794 = 0 

ygg = 0.2057 6131 6872 4279 8353 9094 6502 0576 * 10 3 

y<jg = 0. 1071 4285 7142 8571 4285 7142 8571 4286 

7 97 = 0.3009 9941 2110 5232 2163 4332 7454 4386 • 10 _1 

798 = 0.6666 6666 6666 6666 6666 6666 6666 6667 * 10" 1 

y i00 = 0.4594 4814 6336 7656 7832 1273 3592 1819 • 10" 2 
7ioi = 0 

7102= 0 
7103 = 0 

7ioj =-0. 1740 6947 8190 1750 4798 1220 0310 8640 * 10 -i 

7ioe =-0. 1358 5013 7797 6751 8478 4965 3891 1343 • 10 -3 

yioe = 0. 5624 349 0 8687 4288 9606 4146 7180 2943 * 10“ 3 

7 i07 =-0.4532 9040 6942 0996 9343 6621 7872 6371 * 10 -2 
y 108 = 0. 2139 7111 2333 0803 7669 1281 7128 6734 • 10 -1 

y 109 = 0. 1872 5071 1050 2209 3099 0486 4784 3256 • 10“ 2 

y 110 = 0.2618 1558 1123 7916 3028 4971 0405 9137 • 10 -1 
7m = 0 

7112 = 0 

7113 = 0 

7u4 =-0. 1348 5552 1809 8249 5201 8779 9689 1360 * 10" 1 

y 115 = 0.6168 2138 1305 5131 8451 1711 8731 7478 • 10 -3 

y 116 = 0. 1094 1595 9245 9721 8867 0006 8800 0631 * 10 +1 

7 m = 0.8450 8115 2929 3391 0997 5855 7562 8239 • 10 -2 

7iis = 0. 1777 6087 9419 6241 2575 1394 6991 7887 

7 119 = 0. 6623 0832 0078 1351 3201 1614 2625 9057 • 10 -1 

71110=- 0.9662 6577 5418 7842 5350 4551 0321 1566 
7i2o= 0.1880 6475 4937 1829 2072 2015 3722 2113 * 10 -1 
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TABLE 3. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 8(9)-13-x 

(Concluded) 


7121 

= 

0 


7122 

= 

0 


7123 

= 

0 


7124 

= 

0 


7125 

= 

0 


7126 

= 

0. 

25 

7127 

= 

0. 

3047 

7128 

= 

0. 

1905 

7129 

= 

0. 

8544 

7 121 0 

= - 

-0. 

2344 

71211 

= 

0. 

1967 

7130 = 

C 0 = 

0. 

3214 

7131 = 

:C 1 = 

0 


7132 = 

=C 2 = 

0 


7133 = 

:C 3 = 

0 


7134 = 

:C 4 = 

0 


7 135 = 

IC 5 = 

0 


7136 = 

= c 6 = 

0 


7137 = 

= c 7 = 

0. 

5079 

CO 

CO 

11 

=c 8 = 

0. 

1577 

7139 = 

- c 9 = 

0. 

,7889 

7l310F c 10= 

0. 

, 1600 

71311= 

=°n= 

0. 

,2032 


1545 4235 3444 
2696 2749 1849 
3054 8874 9933 
9246 2505 6950 
6586 5030 1368 
2857 1428 5714 


3650 7936 5079 
9220 7792 2077 
6103 8961 0389 
4659 7153 8134 
8583 2213 6693 


9597 5497 3472 
7140 7193 6294 
1943 7763 8822 
7305 3300 7835 
5512 2865 4566 
2857 1428 5714 


3650 7936 5079 
9220 7792 2077 
6103 8961 0389 
3587 9856 8795 
9300 5183 0086 


3946 

1585 

0833 

2450 

2269 

2857 


3651 

9221 

6104 

0545 

9726 


10 +1 

10' 1 

10"' 1 
10 +1 

10 -1 
i n- 1 


io _1 

10 _1 

10 _1 


74312 =° 12 = 0 


°0 

= 

0. 3214 

2857 

1428 

5714 

2857 

1428 

5714 

2857 * 

°1 

= 

0 








C 2 

= 

0 








C 3 

= 

0 








C 4 

= 

0 








C S 

= 

0 








C 6 

= 

0 








c 7 

= 

0. 1015 

8730 

1587 

3015 

8730 

1587 

3015 

8730 

C 8 


0. 2366 

8831 

1688 

3116 

8831 

1688 

3116 

8831 

C 9 

— 

0. 2366 

8831 

1688 

3116 

8831 

1688 

3116 

8831 

c 10 

— 

0. 1803 

7518 

0375 

1803 

7518 

0375 

1803 

7518 

Cll 

— 

0. 1803 

7518 

0375 

1803 

7518 

0375 

1803 

7518 

Cl2 

= 

0. 3214 

2857 

1428 

5714 

2857 

1428 

5714 

2857 • 
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TABLE 6. COEFFICIENTS (IN FRACTION FORM) FOR RKN 7(8)-10-x 



^ = I^( 531 - 73 ^) - V 61 = 0 , y 62 = 0 , y 63 = _ Ti -Z^( 10 5 - 73^2T) 


y « = i 5 ^( 399+ 73 ^) ’ = u^C 21 - 5 ^) 

= ifliss ( 789 - 89 ^I) > y n - 0 . Y u = o , y 7J = - g~J (mb - «'/«) 

y » = - ( 77 + 103 ’ y ” = who ( 147 - 29 = ITois ( 325 - 51 ^ 21 ) 

>80 = 29l00 ( 5?9+ 235 ' /TI ) ’ y « “ 0 • >« = 0 • y »3 = ~ (315 - 191^21) . Vw = ~ Ti^ ( 161 - 15 

>® = -so ( 21 - 5 ^i) • = irki 395451 ^ 11 ) • y * = m( 43 + 9sfYl ) 

14 1 8 64 

790 = 20 - 49 Y&3 ’ 791=0 * y 82 = 0 , T93 = T» . y*i = - g Y93 . Y* = 0 * 7&e ~ 45 + 441 * 793 

>« - 5^( 7+ ^ 2i ) - i£( 21+5 ^K ■ - 5^( 7 -' /Tl ) - iki 21 - 5 ^)^ 

rioo = <=0 = L , y 101 = C, = 0 , r ,02 = Cj -- 0 , 7,03 = c 3 = 0 - I'm = c 4 = 0 " >105 = C S = 0 

Yik = c 8 = • Tin = c 7 = 553 (7 + -sT2l) , y,08 = <=8 = TJJJ (7 -V2I) , Tw = c 9 = 0 

1 . . . . 16 . 49 . _ 49 • _ 1 

C ° = 20 ’ c < = 0 ’ <= 2 =0, c 3 = 0 , c 4 = 0 , c 5 = 0 , c 6 = - , c, = — , <=8 - ■ c 9 " 20 

Truncation Error in x: TE^. = — (f 9 - f t0 ) h 
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TABLE 7. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 7(8)-10-x 


a, = 0.2380 

9523 

8095 

2380 

9523 

8095 

2380 

9524 



a 2 = 0.4761 

9047 

6190 

4761 

9047 

6190 

4761 

9048 



a 3 = 0. 1428 

5714 

2857 

1428 

5714 

2857 

1428 

5714 



a A = 0. 7142 

8571 

4285 

7142 

8571 

4285 

7142 

8571 



= 0. 8273 

2683 

5353 

9885 

7189 

9146 

2281 

2343 



a 6 = 0. 5 










a 1 = 0. 1726 

7316 

4646 

0114 

2810 

0853 

7718 

7657 



a 8 = 0.8273 

2683 

5353 

9885 

7189 

9146 

2281 

2343 



« 9=1 










«io= 1- 










•y 10 = 0.2834 

4671 

2018 

1405 

8956 

9160 

9977 

3243 

• 

10" 1 

•y 20 = 0 . 3 7 79 

2894 

9357 

5207 

8609 

2214 

6636 

4324 

* 

10" 1 

y 2 i = 0.7558 

5789 

8715 

0415 

7218 

4429 

3272 

8647 

• 

10" 1 

730= 0.7448 

9795 

9183 

6734 

6938 

7755 

1020 

4082 

* 

10~ 2 

7 31 = 0. 3469 

3877 

5510 

2040 

8163 

2653 

0612 

2449 

• 

10" 2 

Y 2 2 = ”0. 7142 

8571 

4285 

7142 

8571 

4285 

7142 

8571 

• 

10“ 3 

740 = 0.2 125 

8503 

4013 

6054 

4217 

6870 

7482 

9932 

* 

10" 1 

741= 0 










7 42 = 0.8199 

7084 

5481 

0495 

6268 

2215 

7434 

4023 

• 

10 _1 

7 43 = 0. 1518 

4645 

2866 

8610 

3012 

6336 

2487 

8523 



7 50 = 0 . 23 7 0 

3048 

2317 

0896 

0011 

4506 

5345 

5765 

• 

10 _1 

751= 0 










7 52 = 0.1178 

3967 

9152 

0302 

1396 

2859 

5090 

6980 



II 

O 

H* 

00 

CO 

00 

9983 

9035 

7159 

7291 

5722 

6756 

6178 



7 54 = 0. 16 79 

2279 

8289 

6771 

0488 

1117 

0444 

5794 

* 

10" 1 

7 e o= 0. 1023 

2915 

3264 

6998 

3308 

2850 

2889 

8847 

• 

10" 1 

7 6 i= 0 










7 62 = 0 










7 63 = 0. 1046 

0261 

5893 

3876 

4605 

2505 

9971 

3509 



7 64 = 0.2228 

6007 

7262 

2584 

3070 

1670 

6647 

5673 

• 

10” 1 

7 g5 =-0. 1212 

1538 

9460 

8347 

2430 

9580 

9250 

9611 

• 

10 _1 

7 70 = 0. 12 9 6 

4311 

6717 

3232 

1068 

4919 

6597 

8170 

• 

10 _1 

7n= 0 










7 72 = 0 










y 73 =-0. 3225 

7396 

6497 

0323 

7961 

9085 

6869 

8362 


10 -2 

7 74 =- 0. 2178 

5924 

4674 

7823 

4947 

5622 

5638 

7063 


10 -1 

7 7S = 0. 1039 

5065 

9297 

9865 

7455 

5729 

1428 

3704 


10" 1 

7 76 = 0. 1656 

0297 

4253 

5186 

5698 

6865 

4750 

4739 


10 _1 

7 8 o= 0.5632 

3309 

1263 

4770 

0732 

9316 

6974 

5785 


10" 1 

7 8 i= 0 
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TABLE 7. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 7(8)-10-x 

(Concluded) 


r 8 2 = 

0 









^83 

■0. 1333 

9808 

5175 

3727 

2410 

9123 

0985 

7192 


T84 =- 

0. 7322 

3305 

2187 

7968 

2461 

2534 

0642 

3873 ' 

10 -2 

Tss = 

0. 1171 

1500 

8659 

9510 

2242 

4912 

4062 

4327 * 

10 -i 

rae = 

0. 1140 

5194 

7472 

6073 

1797 

4782 

8402 

5036 


V87 = 

0.3008 

6850 

4480 

7234 

2878 

3187 

2312 

2697 


11 

<=> 

£ 

0.4183 

6734 

6938 

7755 

1020 

4081 

6326 

5306 • 

10" 1 

Y#i = 

0 









Y 92 = 

0 









II 

£ 

0. 1 









794 =- 

■ 0. 1111 

1111 

1111 

1111 

1111 

1111 

1111 

1111 • 

10' 1 

795 

0 









796 = 

0. 1922 

9024 

9433 

1065 

7596 

3718 

8208 

6168 


797 = 

0. 1256 

4106 

1472 

7145 

7601 

7125 

3258 

9020 


798 = 

0.5134 

3065 

5114 

1240 

8109 

8588 

0109 

3928 • 

10 _1 

7ioo" c (r 

0. 5 







* 

10- 1 

7ioi =c i = 

0 









7l02 :=c 2 = 

0 









7l03 =c 3^ 

0 









7l04~ c 4 = 

0 









7 105 =c 5 = 

0 









11 

CD 

O 

11 

£ 

0. 1777 

7777 

7777 

7777 

7777 

7777 

7777 

7778 


7l07 =c 7 = 

0. 2252 

1674 

9624 

1413 

3346 

1434 

2509 

8916 

10 _1 

7l08~ c 8 = 

0. 4700 

5472 

5980 

8088 

8760 

7879 

7123 

3067 • 

7l09“ c 9 = 

0 









Co 

0.5 







• 

10- 1 

Cj 

0 









c 2 

0 









C 3 

0 









c 4 

0 









C S 

0 









^6 

0. 3555 

5555 

5555 

5555 

5555 

5555 

5555 

5556 


c 7 

0. 2722 

2222 

2222 

2222 

2222 

2222 

2222 

2222 


C 8 

0. 2722 

2222 

2222 

2222 

2222 

2222 

2222 

2222 


C9 

0. 5 







• 

10" 1 













TABLE 9. COEFFICIENTS (IN FRACTION FORM) FOR RKN 6(7)-8-x 


1 

= 5 ’ 

2 

«2 " 5 . 

1 

a* = — , 

3 6 

1 

“4 = 2 * 

tt5= i(7-V2l) , 

«« = n ( 7+ ^r) > 

«7 = 1 . 

1 «8 = 1 






1 



2 

4 





r 2 o = 

75 * 

Y 21 ?5 






277 


95 


35 


Yso = 

31 104 

» Y 31 = 

15 552 ’ 

Y 32 " 

31 104 



5 



25 


9 

Y 40 = 

192 * 

0 

II 

<2 

£ 

II 

1344 

» Y 43 = 

112 


• y5i = 0 ■ = -wki 77 - 16 ^ 1 ) • y^ = ^k( 133 - 27 ^) 

y* = THt (441 - 95 *^ 1 ) 

= irfeo ( 781+ 103 ^T) ’ ** = 0 • *« = - 279^8 (^9 + 599^) 

ro =- 5So ( 2389 + 513 " /Fl ) * y « = iiio ( 315 + 127 ^ = Jb ( 225 + 49 ' V * i ) 


Y 70 = 

20 ” 

_L_ 

54 

Y 73 » 

Yn ” 

0 , 

Y 72 ” 0 » 

Y73 

" Y 73 ♦ 

ira 

^|oo 

+ 

-IS 

II 


Yts = 

360 1 

I 7 + 

V 21 ) 

--Z-I 

324 ' 

^23 + 5 

'JH.)y n , 

y 76 

= — ( 
360 V 

7 -^ 2i )-H4 ( 23 

- 5 V 2 l)y 73 

Yso = 

c 0 = 

J_ 

20 

. Ybi = c t 

= 0 

> Y82 " C 2 

= 0 

» Y83 

II 

£ 

0 

II 

C5 

O 

II 

£ 

11 

£|°° 

Ysb = 

c 5 = 

— (7 + 

360 ' 

^2l) 

. Yge 

7 

° e 360 

( 7 - 

V 21 ), 

3* 

II 

0 

II 

0 



1 







16 

49 

49 

Co = 

20 * 

Cl 

= 0 

, c 2 

= 0 , 

. c 3 = 0 , 

C4 

45 ’ 

J? 

II 

II 

0 

o» 

180 * C? 


Truncation Error in x: TE^. = (fj - f 8 ) h 
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TABLE 10. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 6(7)-8-x 


04= 0. 2 










a 2 = 0.4 










o^3 - 0, 1666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 



« 4 = 0.5 










a 5 = 0.1726 

7316 

4646 

0114 

2810 

0853 

7718 

7657 



a 6 = 0.8273 

2683 

5353 

9885 

7189 

9146 

2281 

2343 



1 










a 8 = 1 










7i 0 ” 0. 2 








• 

10" 1 

72o = 0* 2666 

6666 

6666 

6666 

6666 

6666 

6666 

6667 


10" 1 

y 2 i= 0.5333 

3333 

3333 

3333 

3333 

3333 

3333 

3333 


10“ 1 

r 3 0 " 0* 8905 

6069 

9588 

4773 

6625 

5144 

0329 

2181 


10“ 2 

y 31 = 0.6108 

5390 

9465 

0205 

7613 

1687 

2427 

9835 


10~ 2 

732—0. 1125 

2572 

0164 

6090 

5349 

7942 

3868 

3128 


10~ 2 

740= 0. 2604 

1666 

6666 

6666 

6666 

6666 

6666 

6667 


10 _1 

T 4 i= 0 










y 4 2= 0. 1860 

1190 

4761 

9047 

6190 

4761 

9047 

6190 


10 -1 

7 43 = 0. 8035 

7142 

8571 

4285 

7142 

8571 

4285 

7143 

• 

10" 1 

y 50 = 0.8038 

6592 

6268 

7269 

1853 

8866 

9589 

7376 

• 

10" 2 

Tsi" 0 










752—0.3192 

0630 

9932 

1949 

0963 

7585 

4592 

1420 

• 

10 -2 

753 = 0.8687 

4329 

5769 

7925 

6975 

3275 

0096 

2197 

• 

10 -2 

7 54 = 0. 1373 

9817 

7337 

1039 

6924 

5274 

9415 

8987 

• 

10 -2 

760= 0. 1065 

4806 

9437 

1132 

2454 

7497 

3521 

2194 



r S i= 0 










y 62 =-0. 3674 

8509 

6867 

1886 

7582 

2671 

9499 

4454 



7 63 — 0. 6218 

4769 

6554 

0978 

6166 

7886 

4999 

0441 

• 

10 +1 

764= 0.3050 

9765 

7571 

2216 

6014 

8531 

2903 

4131 



y 6 5= 0. 6516 

5311 

8164 

8255 

2651 

8365 

2849 

1585 

• 

10 +1 

y 7 o= 0.4981 

4814 

8148 

1481 

4814 

8148 

1481 

4815 

• 

10“ 1 

771= 0 










772= 0 










773= 0. 1 








• 

10“ 1 

774= 0. 1779 

0123 

4567 

9012 

3456 

7901 

2345 

6790 



775 = 0 . 2 152 

9730 

0570 

9482 

9641 

7280 

4963 

0457 



776= 0.4698 

6650 

0463 

3565 

4200 

0034 

5431 

2711 

• 

10 -1 
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TABLE 10. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 6(7)-8-x 

(Concluded) 


Tso-Co = 

T8i =c r 

T 82 =c 2 = 

T83 =c 3 = 

T84 =c 4" 

?85 =c 5" 

'>' 86 == c 6 = 

y%f zc f 


0.1777 7777 7777 7777 7777 7777 7777 7778 
0.2252 1674 9624 1413 3346 1434 2509 8916 
0.4700 5472 5980 8088 8760 7879 7123 3067 
0 


-0. 3555 5555 5555 5555 5555 5555 5555 5556 
=0. 2722 2222 2222 2222 2222 2222 2222 2222 
=0.2722 2222 2222 2222 2222 2222 2222 2222 
=0.5 


TABLE 11. PATTERN FOR RKN 5(6)-6-x 


2 3 4 5 


6 















TABLE 12. COEFFICIENTS (IN FRACTION FORM) FOR RKN 5(6)-6-x 
“ 1 = 1 . “2 = JZ ( 5 + ' /T ) • “8 = 13 (5 - . 

“* = To ( 5 + ’ “ 5 = 1 ' “e = 1 


1 


7io = 

2 







T 20 = 

1 

300 

(35 + 

11 • \l~tj 

, 7 2 i 


• nTs) 

7 3 o = 

1 

600 

(25 - 

3 \ ft) , 

7 3 i = 

_ _L_ sfT 
300 

r 32 = ^0 («-5VI) 

T 40 = 

1 

600 

(25 + 

3 n/T) , 

7u = 

1 

300 

n/T , 

742 = 0 

743 = 

1 

120 

(13 + 

5 n/T) 





750 = 

_1_ 
12 " 

• 751 

* 751 = 

7 5 i , 

7 52 

= 0 , 


753 c 

-±-( 
24 V 

5 + Nf5 ) + ‘JT- 

751 » 

754 

-£(• 

- i n/ _ 5^ - nT 5 * 7si 

7 6 o z 

= c 0 = 

1 

12 

. 7«i = 

c 4 = 0 

> 

762 = c 2 

= 0 , 

7 S 3 

= C 3 = 

24 ' 

^5 + \/~5) 

» 7 64 

= c 4 

-£(• 

- nTs) , 7 65 = c 5 = 0 


1 # .5.5,1 

c 0 ~ ^2 9 C 1 — ® > c 2 “ ® > ^3 “ ±2 9 ^4 — Y2 9 ~ 12 


Truncation Error in x: TE. = -pr- (f 5 - f 6 ) h 

X xz 
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TABLE 13. COEFFICIENTS (IN DECIMAL FORM) FOR RKN 5(6)-6-x 


«l = 1 

Q! 2 = 0.7236 0679 7749 9789 6964 0917 3668 7313 

<*3 = 0.2763 9320 2250 0210 3035 9082 6331 2687 

a A = 0.7236 0679 7749 9789 6964 0917 3668 7313 

«6 = 1 

«e = 1 

Tio = 0. 5 

y 20 = 0. 1986 5582 5841 6589 5553 5003 0345 2015 

y 21 = 0.6314 7573 0333 3052 9285 4556 4891 6417 •- 10 -1 

y 30 = 0.3048 6326 7791 6771 8184 6207 9832 3010 • 10 _i 

y 31 =-0.7453 5599 2499 9298 9880 3057 8895 7709 • 10"2 

y 32 = 0. 1516 3834 2708 4209 5982 9510 9713 6197 • 10 _1 

y 40 = 0.5284 7 006 554 1 656 1 5148 7125 35 01 0323 * 10 _1 

y 41 = 0.7453 5599 2499 9298 9880 3057 8895 7709 * 10" 2 

T42 = 0 

y 43 = 0.2015 0283 2395 8245 7068 3715 5695 3047 

y 50 =0.8 • 10-1 

y 51 = 0.3333 3333 3333 3333 3333 3333 3333 3333 • 10 -2 

T52 = 0 

y R3 = 0.3089 5639 2320 8238 6967 1746 1484 2624 

y 54 = 0. 1077 1027 4345 8427 9699 4920 5182 4043 

y 60 =c 0 = 0.8333 3333 3333 3333 3333 3333 3333 3333 • 10 _1 
y 6i =Ci= 0 
7 / S2 =c 2 = 0 

y 63 =c 3 = 0.3015 0283 2395 8245 7068 3715 5695 3047 
yg 4 =e 4 = 0. 1151 6383 4270 8420 9598 2951 0971 3620 

'y65 =c 5~ 0 

c 0 = 0.8333 3333 3333 3333 3333 3333 3333 3333 • 10 -i 

<$! =0 
c 2 = 0 

6 3 = 0.4166 6666 6666 6666 6666 6666 6666 6667 

c 4 = 0.4166 6666 6666 6666 6666 6666 6666 6667 

c s = 0.8333 3333 3333 3333 3333 3333 3333 3333 * 10 -1 
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TABLE 14. APPLICATION OF THE VARIOUS FORMULAS TO ORBIT NO 
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TABLE 15. APPLICATION OF THE VARIOUS FORMULAS TO ORBIT NO, 
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Figure 1. Central orbits of Section V. 



